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Abstract 

I present a currents algebra for the two-sites Bose-Hubbard model, generalize the Heisen¬ 
berg equation of motion to write the second time derivative of the currents operators and 
use it to get the quantum dynamics of the currents. For different choices of the Hamiltonian 
parameters I get different currents dynamics and determine the period of the oscillations in 
function of the parameters. 


1 Introduction 


The early experimental realization of a two-wells Bose-Einstein condensate (BEC) was made 
only two years after the experimental verihcation of the BEC [IHIl] to study the interference 
between two freely expanding condensates [TSlfTG] . and their results had direct implications 
in the study of the atom laser and the Josephson effect [T^HH] for BEC. Some models were 
used to study some behaviors of these systems as for example the quantum phase transitions, 
the classical analysis and the quantum dynamics [19H23]. I am considering here the two- 
sites Bose-Hubbard model, also known as the Canonical Josephson Hamiltonian [8], that 
has been an useful model in understanding tunneling phenomena using two BECs |2lH30]. 
This model is integrable in the sense that it can be solved by the quantum inverse scattering 
method (QISM) [MfITT] and it has been discussed in different ways using this method [33fH0] . 
In this context this model is a particular case of the bosonic multi-state model studied 
in m The experimental quantum dynamics and the classical analysis of this model was 
performed by |48] - [45] . In this letter I will discuss the currents algebra for the two-sites Bose- 
Hubbard model. Currents algebra was introduced by M. Gell-Mann in high energy physics 
to study partially conserved axial vector current in the beta decay |16]. I generalize the 
Heisenberg equation of motion to write the second time derivative of any operator and use 
it to study the quantum dynamics of the currents. This method can be applied to many 
systems that present microscopic tunneling phenomenon to get some characteristic energies 
of the systems in function of the period of the oscillation and that is also important to 
quantum metrology |47H53]. The model is described by the Hamiltonian 

H = f (iVi - N^f - - N^) - ^{a\a2 + a^ai), (1.1) 

o Z Z 

where a\, al, denote the single-particle creation boson operators in the two wells and Ni = 
a|ai,iV 2 = 0 , 20 , 2 , are the corresponding number of particles boson operators. These bosons 
operators satishes the canonical commutation relations [aj,aj] = SijI, [ai,aj] = [«!,«]] = 0 
and [Ni,aj] = —SijOj, [iVj,a]] = HhijOj, where I is the identity operator. The coupling K 
provides the strength of the s-wave scattering interaction between the bosons, A/i is the 
external potential and Sj is the amplitude of tunneling. 
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2 Symmetries 


The Hamiltonian (II.Ih is invariant under the Z 2 mirror transformation hj aj 

and under the global U{1) gauge transformation cij —)■ e^°‘aj, where a is an arbitrary c- 
number and aj —)■ j = 1,2. For a = n we get again the Z 2 symmetry. The global 

f/(l) gauge invariance is associated with the conservation of the total number of atoms 
= A^i + A '2 and the Z 2 symmetry is associated with the parity of the wave function by 
the relation P |\ 1 /) = (—l)^|\h), with 


^ vn! 


V ... 

where P is the parity operator and [H, P] = 0. 


| 0 , 0 ), 


( 2 . 2 ) 


There is also the permutation symmetry of the atoms of the two wells if we have A/i = 0 
and when we turn on Ap we break the symmetry. The wave function 02 . 21 ) is symmetric 
under this permutation 


N 


vrn^Y. 


N—n,n 




( 4 )’ 


n=0 


^/w 


n] 


y/nJ. 


|0,0) = |vl/), 


(2.3) 


where V is the permutation operator and [P, P] = 0 if Afi = 0 [33]. In the antisymmetric 
case A/i 7 ^ 0 we can change the bias of one well. In this case it is called a tilted two-wells 
potential |2HIIS]. In the Fig. ([T|) we represent the two BECs by a two-wells potential for 
the case A/r 7 ^ 0. We get the two-sites Bose-Hubbard model when we consider each BEC as 
a site. 


3 Currents Algebra 


The total particles number boson operator, N = Ni + N 2 , is a conserved quantity and it is 
a commutable compatible operator with the particles number bosons operators in each well, 
[iV, iVi] = [iV, iV2] = [Pi,iV2] = 0 . The number of particles bosons operators in each well 
don’t commute with the Hamiltonian and their time evolution is dictated by the Josephson 
tunneling current operator, = ^(a|a 2 — a\ai), in coherent opposite phases because of the 
conservancy of N. We get the following equations for the time evolution 
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Figure 1: (Color online) A two-wells potential representation of the two-sites Bose-Hubbard 
model for the case A/i 7 ^ 0 and barrier height Vq. We are considering one condensate, in 
light blue, for each well. We show one atom, in dark blue, tunneling from left to right (red 
arrow). 


iVi(t) = W(0)-^^ i(r)dr, (3.4) 

Mt) = Mo) + ^ I ( 3 - 5 ) 

Here is worth to note that the two BECs are entangled by the tunneling of the particles 
and we can study the quantum phase transition of the system using tools of the quantum 
information isniEi]. 

The tunneling current together with the imbalance current X = |(iVi — W) and the 
coherent correlation tunneling current operator T = ^{a\a 2 + aldi), generate together the 
currents algebra [T, jf] = +ii, [T, X] = and [iX,X] = +iT'. With the identihcation 
X 3 . = hT, Ly = hjf, and = hX we can write this currents algebra in the standard compact 
way of the momentum angular algebra [X^., Li\ = ihekimLm, where Skim is the antisymmetric 
Levi-Civita tensor with k,l,m = x, y, z and Sxyz = +1- We have two Casimir operators for 
that currents algebra. One of them is the total number of particles, €i = N, related to the 
U{1) symmetry and the another one is related to the momentum angular algebra and the 
0(3) symmetry, C 2 = -l-X^ -l- iX^. It is direct to show that €2 is just a function of £1 
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and that the Casimir surface is spherical with radius (£ 2 )- 


4 Currents Quantum Dynamics 


We can rewrite the Hamiltonian fll.ip using the currents operators 

H = - Sjf. 


(4.7) 


The quantum dynamic of the currents are determined by the currents algebra, their 
commutation relations with the Hamiltonian and the parameters. If the Hamiltonian is 
not explicitly time-dependent it is not time-dependent, ^ ^ = 0 ) cind the system is 

closed (conservative). It is also important to note that the Hamiltonian is the same in the 
Heisenberg and Schrodinger pictures, Hh = Hs- Using these facts we can write the second 
time derivative of any operator O in the Heisenberg picture as 

(4.8) 


or as 

d^O i ' dO 
dt‘^ h ’ dt 


(4.9) 


It is direct to generalize the Eqs. fl4.8p and fl4.9p for higher time derivatives. We can get 
similar equation in the another pictures. The pictures preserve the commutation relations 
between the operators in the sense that if we have [A 5 , Bs] = Cs in the Schrodinger picture 
we get the same relation in the Heisenberg picture, [A^, B^] = Ch, and in the interaction 
picture, [Aj, Bj] = Cj. The same is true for the anticommutators, and so the pictures 
preserve the algebra. We can see from Eqs. fl4.7p and fl4.8p that the Casimir operators Cfi 
and £2 are also conserved quantities, [H, £ 1 ] = [H, £ 2 ] =0. Using the Eq. fl4.8p or fl4.9p we 
found the following equations for the quantum dynamics of the three currents 


d^X 

-h —X 

dB 


S.jK 


h? 


XT 


.SjK - 

—— A'T — 7 


(4.10) 


dB 
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(A/i)^ + £j + 


~T 


J 




K8j , KA^Ji^ i KSj^ 


(4.11) 
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J. (4.12) 


We can see from Eqs. fl4.injl . fl4.1ip and fl4.12p that the cnrrents are conpled on the right 
hand side of these eqnations. Different choices of the ratio between the parameters of the 
Hamiltonian gives us different dynamics for the currents. In the Rabi regime, K/£j N~‘^ 
[8l|26l|38]. Consequently, in the extreme Rabi regime we can neglect K and consider the 
no interaction limit K ^ 0. Considering the symmetric case, A/i = 0, the current T is a 
conserved quantity, [H, T] = 0, but this don’t means that we don’t have tunneling. We can 
see from Eqs. fl3.4p and fl3.5p that the quantum dynamic of Ni, N 2 , and X only depend of the 
current J' and the amplitude of tunneling Sj. Here is worth to note that in Sj is included 
the kinetic energy of the atoms. The current dynamics for these currents are the dynamic 
of the simple harmonic oscillator (SHO) 


= 0, (4.13) 

= 0, (4.14) 

where ux = ojj = ujsj = ^ is the natural angular frequency of the SHO. The period of the 
oscillations is T = -A_ we get the following relation £jT = h between the energy of the 

^£j 

amplitude of tunneling and the period. In the no interaction limit is expected a period of 
T = 500 ms instead the period of T = 40.1 ms for the interacting nonlinear regime as in the 
experiment [43] and in the generalized model [55] . We have gotten the value 8j = 12.5h J 
for the amplitude of tunneling. In the Fig. ([2|) we show the solutions for the Eqs. fl4.13|) 
and fl4.14p . The currents are uncorrelated now and there is no interference between them. 
We have Rabi dynamics for the currents X and JT and self-trapping for the current T. 

Breaking the symmetry, A/i 7 ^ 0, to consider the antisymmetric case the currents dynamics 
are 


d ‘^1 


dt'^ 

d?J 

dt‘^ 


+ 




d^X f£j\\ 


(4.15) 
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Figure 2: Quantum dynamics of the currents for ojx = oJj = 12.5 rad-Hz. The initial 
condition for the current X{t) (full line) is X(0) = 1.0. The initial condition for the current 
J'{t) (dashed line) is 17(0) = —1.0. The initial condition for the current T{t) (dot-dashed 
line) is T(0) = 0.4. The initial conditions for the hrst derivative of all currents is zero. 


d^f 


(4.16) 


The Eq. (I4.17|l describes a SHO with natural angular frequency ojj = ^, period of 
the oscillations T = , —= and relation ET = h, with E = hwj, between the period and 

the energy. The Eqs. fl4.15p and fl4.16p are a system of two second order linear differential 
equations. If we diagonalize the matrix of the coefficients we get the same frequency uj. In 
the Fig. ([3]) we show the numeric solution for the Eqs. fl4.15l) . fl4.16l) and fl4.17l) . We choose 
the same period, T = 500 ms to get the values £j = 12.5h J for the amplitude of tunneling 
and A/i = 0.5h J for the external potential. We have Rabi dynamics for the currents J7, 
Josephson dynamics for the current X and self-trapping for the current T. The currents X 
and T are correlated and there is interference between them. 


6 
































t{ms) 


Figure 3: Quantum dynamics of the currents for ujgj = 12.5 rad-Hz and = 0.5 rad-Hz. 
The initial condition for the current X{t) (full line) is X(0) = —1.0. The initial condition 
for the current J{t) (dashed line) is iT(0) = 1.0. The initial condition for the current T(t) 
(dot-dashed line) is T(0) = —0.9. The initial conditions for the first derivative of all currents 
is zero. 


5 Summary 

In summary, I have showed that a currents algebra appears when we calculate the quantum 
dynamics of the number bosons operators of each well. I have generalized the Heisenberg 
equation of motion to write the second time derivative of any operator. Then I have cal¬ 
culated the quantum dynamics of these currents and have showed that different dynamics 
appear when we consider different choices of the parameters of the Hamiltonian. For specihc 
choices of the parameters some of the currents are uncorrelated and there is no interference 
between them. 
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